INVERSE PROBLEM IN CYLINDRICAL ELECTRICAL NETWORKS 
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Abstract. In this paper we study the inverse Dirichlet-to-Neumann problem for certain 
cylindrical electrical networks. We define and study a birational transformation acting on 
cylindrical electrical networks called the electrical i?-matrix. We use this transformation 
to formulate a general conjectural solution to this inverse problem on the cylinder. This 
conjecture extends work of Curtis, Ingerman, and Morrow |CIMj . and of de Verdiere, 
Gitler, and Vertigan [dVGVj for circular planar electrical networks. We show that our 
conjectural solution holds for certain "purely cylindrical" networks. Here we apply the 
grove combinatorics introduced by Kenyon and Wilson |KWj . 



1. Introduction 

In this paper we consider the simplest of electrical networks, namely those that consist 
of only resistors. The electrical properties of such a network N are completely described 
by the response matrix L(N), which computes the current that flows through the network 
when certain voltages are fixed at the boundary vertices of N. 

De Verdiere-Gitler- Vertigan [dVGV] and Curtis-Ingerman-Morrow [CIM] studied the 
inverse (Dirichlet-to-Neumann) problem for circular planar electrical networks. Specif- 
ically, they considered networks embedded in a disk without crossings, with boundary 
vertices located on the boundary of the disk. The following theorem summarizes their 
results. 

Theorem. 

(1) Any circular planar electrical network is electrically equivalent to some critical 
network, which is characterized by its medial graph being lenseless (see |dVGV 
Theoreme 2]) 



(2) Any two circular planar electrical networks having the same response matrix can be 
connected by simple local transformations: series-parallel, loop removal, pendant 
removal, and star-triangle transformations discussed in Section [Ol Furthermore, 
if both networks are critical, then only star-triangle transformations are required 
(see jdVGVl Theoreme 4] or jCLMl Theorem I]). 

(3) The edge conductances of a critical circular planar electrical network can be re- 
covered uniquely from the response matrix (see [CIM] Theorem 2] or [dVGV, 
Theoreme 3]). 

(4) The response matrices realizable by circular planar networks are the ones having 
all circular minors nonnegative (see \CIM\ Theorem 3] ). 

(5) The space Y of response matrices of circular planar networks has a stratification by 
cells Y = L\Ci where each C\ ~ R> can ^ e obtained as the set of response matrices 
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for a fixed critical network with varying edge weights (see [dVGVt Theoreme 3 and 
5] or [ClMl Theorem 4]j. 

It is an open problem to extend these results to electrical networks embedded in surfaces 
with more complicated topology. In this paper we make progress towards understanding 
the inverse problem for networks embedded in a cylinder. Our first main result is to 
construct a birational transformation we call the electrical R-matrix. This transformation 
acts on the edge weights of a local portion of an electrical network embedded into the 
cylinder, preserving all the electrical properties of the network (Corollary I5.4p . Further- 
more, this electrical .R-matrix satisfies the Yang-Baxter relation (Theorem 14 . 5 j) . and is a 
close analogue of the "geometric .R-matrices" of affine crystals, to be explained below. 

Using the electrical .R-matrix, we formulate the following general conjecture. A more 
precise version is given as Conjecture 16.31 

Conjecture. 

(1') Any cylindrical electrical network is electrically equivalent to a critical cylindrical 
electrical network. 

(2') Any two cylindrical electrical networks G and G' with the same universal response 
matrices are connected by local electrical equivalences. Furthermore, if G and G' 
are both critical, then only star-triangle transformations, and electrical R-matrix 
transformations are needed. 

(3') // a cylindrical electrical network is critical, then the edge conductances can be 
recovered up to the electrical R-matrix action. 

(5') The space X of universal response matrices of cylindrical electrical networks has 
an infinite stratification by X = UCj where each Ci ~ M> x M> is a semi-closed 
cell that can be obtained as the set of universal response matrices for a fixed critical 
network with varying edge weights. 

The naive analogue of (4) does not hold - see Section 1531 The universal response matrix 
in the Conjecture is the response matrix of the universal cover of the cylindrical network 
G. Roughly speaking, it allows us to not only measure the current flowing through the 
boundary vertices, but also how many times the current has winded around the cylinder. 
It may be possible to formulate this in a more electrically natural way by measuring 
magnetic fields. 

Thus the key difference between the planar and the cylindrical cases is that even for 
a critical network on a cylinder the edge conductances may not be uniquely determined 
from the response matrix. This non-uniqueness comes from the existence of the electrical 
.R-matrix, the action of which preserves both the underlying graph of the network and 
the response matrix, while changing the edge conductances. The action of the electrical 
.R-matrix can on the one hand be thought of as a Galois group, and on the other hand as 
a monodromy group. 

Our second main result is to establish the above Conjecture for a certain class of 
cylindrical networks we denote N(m). These critical networks can be thought of as the 
"purely cylindrical" networks. There is no local configuration for which the star-triangle 
transformation can be applied in N(m), but the electrical .R-matrix generates an action 
of the symmetric group S m . We show in Theorem 16.21 that for the networks N(m) the 
edge conductances are recovered uniquely from the universal response matrix modulo 
this S m action, and in particular the inverse problem has generically ml solutions. One 
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way to formulate our main conjecture is that the networks N(m) exactly encapsulate the 
difference between the planar and cylindrical cases. 

The proof of Theorem 16.21 occupies the technical heart of the paper: we express the 
edge conductances as limits of certain rational functions of the universal response matrix 
(Theorem 15. lip . Here we use crucially the work of Kenyon and Wilson [KW\ IKW2j . 
Kenyon and Wilson study groves in circular planar electrical networks. These are forests 
whose connected components contain specified boundary vertices. Kenyon and Wilson 
connect ratios of grove generating functions with the response matrix of the corresponding 
network. By a careful choice of grove generating functions, we can recover the desired 
edge conductances. 

Recall that a matrix is totally nonnegative if every minor of it is nonnegative. There is a 
mysterious similarity between electrical networks and a different kind of networks arising 
in the theory of totally nonnegative matrices. In |LPlla+l , we presented an approach to 
understanding this similarity via Lie theory. Whereas the theory of total nonnegative is 
intimately related to the class of semisimple Lie groups ( |Lusj ). we suggested in | LPlla+ 
that a different class of "electrical Lie groups" is related to electrical networks. These 
electrical Lie groups are certain deformed versions of the maximal unipotent subgroup of 
a semisimple Lie group. 

The main ideas of the present work are also motivated by this analogy, though our 
philosophy here is more combinatorial in nature. The construction of the electrical R- 
matrix follows the techniques developed in | LP10+ |. There we constructed, using purely 
network-theoretic methods, the geometric (or birational) i?-matrix of a tensor product of 
affine geometric crystals for the symmetric power representations of U'(sl n ). In this paper, 
we use electrical networks instead of the "totally nonnegative" networks of | LP10+| , but 
nevertheless the underlying combinatorics is developed in parallel. We plan to expand on 
this analogy in |L Pllb+| . 

Our formulation of Conjecture 16.31 is also motivated by the analogy with total non- 
negativity. Indeed, analogues of (l')-(5') (and even the missing (4')!) for the totally 
nonnegative part of the rational loop group are established in | LP08+1 |LP10+ |. In par- 
ticular, in | LP08+ | we studied in detail from the totally nonnegative perspective, the 
networks N(m), or more precisely, their medial graphs. Our solution here to the inverse 
problem for the networks N(m) follows the strategy in | LP08+ |, where edge weights are 
recovered by taking limits of ratios of matrix entries; this approach was originally ap- 
plied by Aissen-Schoenberg- Whitney |ASWj to classify totally positive functions. The 
situation we consider here is technically much more demanding, involving rather intricate 
Kenyon- Wilson grove combinatorics. 

Acknowledgements. 
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2. Electrical networks 
For more background on electrical networks, we refer the reader to |CIMt ldVGV[ IKW] . 



2.1. Response matrix. For our purposes, an electrical network is a finite weighted undi- 
rected graph N, where the vertex set is divided into the boundary vertices and the interior 
vertices. The weight w(e) of an edge is to be thought of as the conductance of the corre- 
sponding resistor, and is generally taken to be a positive real number. A O-weighted edge 
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would be the same as having no edge, and an edge with infinite weight would be the same 
as identifying the endpoint vertices. 

We define the Kirchhoff matrix K(N) to be the square matrix with rows and columns 
labeled by the vertices of N as follows: 

K .. = |E e joins i and j W ( e ) for 1 + J 

\ -Ee incident to i W ( e ) for % = j. 

Let M be a square n x n matrix, and Id {1,2,..., n}. Recall that the Schur comple- 
ment M /Mjj is the square n— \I\ matrix defined to be M\^_j^_j — M\ n ^jjMj\Mi ,[ n ]-/, 
where Mj^k denotes the submatrix of M consisting of the rows labeled by J and the 
columns labeled by K. We define the response matrix L(N) to be the square matrix with 
rows and columns labeled by the boundary vertices of N, given by the Schur complement 

L(N) = K/Ki 

where Kj denotes the submatrix of K indexed by the interior vertices. The response 
matrix encodes all the electrical properties of iV that can be measured from the bound- 
ary vertices. Note that our Kirchhoff and response matrices are the negative of those 
commonly used in the literature. 

2.2. Local electrical equivalences of networks. We now discuss the local transfor- 
mations of electrical networks which leave the response matrix invariant. The following 
proposition is well-known and can be found for example in |dVG V] . 

^<^f^><C - ^> a + b .ee ^> a • b '< - ^> ab /( a + b ) t < z 

b 



— > > • — > 

Figure 1. Series-parallel transformations, loop and pendant removal. 

Proposition 2.1. Series-parallel transformations, removing loops, and removing interior 
degree 1 vertices, do not change the response matrix of a network. See Figured 




Figure 2. The Y — A, or star-triangle, transformation. 
The most interesting local transformation is attributed to Kennelly [Ken]. 
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Proposition 2.2 (Y — A transformation). Assume parameters a,b,c,A,B and C are re- 
lated by 

be „ ac a& 
A = , 5 = , C 



a + b + c a + b + c a + b + c 

or equivalently by 

AB + AC + BC , + + + + 

a = A ' h = B ' C= C • 

Then switching a local part of an electrical network between the two options shown in 
Figure [H does not affect the response of the whole network. 

2.3. Tetrahedron relation for electrical networks. 

Proposition 2.3 (Tetrahedron relation). The sequence ofY — A transformations shown 
in Figure [3] returns the conductances of edges to their original values. 









Figure 3. 

Proof. Direct computation. □ 

Remark 2.1. In the terminology of Kashaev, Korepanov, and Sergeev |KKS] . the propo- 
sition states that the Y — A transformation is a solution of type (e) to the functional 
tetrahedron equation. 

3. Kenyon- Wilson L-polynomials 

3.1. Groves. A planar partition of [n] = {1, 2, . . . , n} is a partition r = {r\, t%, . . . ,Tk} 
such that there do not exist i < j < k < £ such that i, k belong to the same part of r, 
and k, i belong to the same part. 

A circular planar electrical network N is an electrical network embedded into a disk, so 
that the intersection of iV with the boundary of the disk is exactly the boundary vertices 
of N. We suppose that the boundary vertices of iV are exactly [n] , and that these vertices 
are arranged in order around the boundary of the disk. 
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A grove in N is a spanning forest where each connected component intersects the 
boundary. A grove T has boundary planar partition r = {ti, t 2 , . . . , r^} if the connected 
components T = Ul^ of T are such that Tj contains the boundary vertices labeled by 7$. 
The weight of the grove is the product of the weights of its edges. Kenyon and Wilson 
[KW] study the probability Pr(r) that a random grove of N is of type r, where the 
probability of a grove is proportional to its weight. 

Let "uncrossing" denote the partition of [n] into singletons. 

Theorem 3.1. [KW, Theorem 1.1, Lemma 4.1] Let G be a finite circular planar electrical 
network. 

(1) The ratio Pr(r) := Pr(r)/Pr(uncrossing) is an integer- coefficient polynomial in 
the Lij(G), homogeneous of degree n — Mparts of r. 

(2) Suppose t be a planar partition with parts of size at most two. Then the polynomial 
of (1) depends only on Lij(G), for i,j 's which are not isolated parts of r. 

3.2. A bound on certain Kenyon- Wilson polynomials. Recall from [KW] that for 
a partition r, we define 

(i) i,=E n l m 

F {i,j}£F 

where the sum is over spanning forests F of the complete graph, for which the trees of F 
span the parts of r, and the product is over the edges of F. 

Let's recall Rule 1 from [KW, p. 5]. If a partition r is non-planar then one can pick 
a < b < c < d so that a, c belong to one part of r, and b, d belong to another part of r. 
Arbitrarily subdivide the part containing a and c into two sets A and C so that a £ A 
and cGC, and similarly obtain B and D. Denote the remaining parts of the partition r 
by "rest". Then the rule is 

AC | RD | rest ->■ A\BCD\rest + £|AOD|rest + C\ABD\iest + D\ABC\rest 

-AB\CD\rest - AD|£C|rest. 

Iterating this rule transforms each non-planar partition r into a linear combination of 
planar ones, and Kenyon and Wilson show that the coefficients do not depend on how 
Rule 1 is applied. Denote by P a%T the coefficient of a planar partition a in the application 
of Rule 1 to t. 

Theorem 3.2 ( |KW1 Theorem 1.2]). We have Pr(<r) = ^ t P <TjT L t . 

Lemma 3.3. Let S be a cyclic interval. Suppose that a planar partition a is such that 
no part p of a contains two elements of S. Let r be a possibly nonplanar partition such 
that a occurs in the expansion of r under repeated application of Rule 1 ( until no more 
applications are possible). Then r does not contain any part p such that \p fl S\ > 2. 

Proof. By |KW[ Theorem 1.2], the result of repeatedly applying Rule 1 does not depend 
on the choices made when applying Rule 1. Let us suppose that r contains a part p such 
that \p fl S\ > 2. In applying Rule 1, if we ever encounter that a, c G p, we will first 
try to choose A and C so that |A Pi 5"| > 2 or \C fl S\ > 2. This would guarantee that 
all the partitions occurring in Rule 1 contain some part p' such that \p' fl S\ > 2. This 
choice is impossible only if a, c G S, which would in turn imply that b or d lies in S. But 
in this case again all the partitions occurring in Rule 1 contain some part p' such that 
\p' H S\ > 2. Thus a cannot occur in the expansion of r. □ 
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Lemma 3.4. Suppose a is a partition such that the non-singleton parts of a contain at 
most K elements. Then there is a constant ck such that when Pr(cx) is expanded as a 
polynomial in the Ly 's the coefficient of any monomial in the 's is less than or equal 
to Ck- 

Proof. Denote the set of boundary vertices by S, and by T C S the elements in the non- 
singleton parts of a. Applying Lemma 13.31 we see that every r, such that a occurs in 
the Rule 1 expansion of r, have parts that contain at most one element from each cyclic 
component of S — T. Since the number of cyclic components is bounded by K, we see that 
there is a bound, depending only on K, on the number of combinatorial types of possible 
r's. Here combinatorial type means the partition one obtains when singletons in S — T 
are removed, and only the relative orders of the remaining elements are remembered. It 
follows that one can find a constant ck so that the coefficient of a in the Rule 1 expansion 
of t is bounded by Ck, for any r. 

Next we note that for the sum of ([I]), the edges of F determine F, which in turn 
determines r. So each monomial in the Ly's occurs in at most one L T . Thus the coefficient 
of any monomial in the Ljj's is bounded by ck- D 

4. Electrical .R-matrix 

The ideas of this section follow closely the calculation of the "whurl transformation" 
in | LP10+ |. In a special case the whurl transformation reduces to the birational, or 
geometric R-matrix of certain affine geometric crystals. This motivates the terminology 
of an electrical R-matrix. 

4.1. The ^-transformation. Fix n > 1 and m > 1. We define an electrical network 
denoted N(m) which is embedded in a cylinder. It has 2n boundary vertices all lying 
on the boundary of the cylinder, with 1,2, ... ,n on the left boundary component and 
1', 2', ...,n' on the right boundary component. There are (m — l)n internal vertices, 
denoted where 1 < j < m — 1 and 1 < i < n. For convenience, the boundary 

vertices 1,2, ... ,n are denoted Adj® and the boundary vertices 1', 2', . . . , n' are denoted 
M^ m \ For each % — 1, 2, . . . , n and j = 0, 1, . . . , m — 1, we have edges from Mj to M^ +1 \ 
and from to M^ l \ Here all lower indices are taken modulo n. 




Figure 4. The cylindrical network N(m). 

Now we focus on the network N(2). We label the edge weights of iV(2) as follows. For 
each i = 1, 2, . . . , n, we have edges with weights <2j from % to M^_ x , weights hi from % to 

M\ ; , weights q+i from to i', and weights di from to (i — 1)'. Here all indices 
are taken modulo n. 



THOMAS LAM AND PAVLO PYLYAVSKYY 

Now define polynomials Ki(a, b, c, d) and Tj(a, b, c, d) as follows: 

n—1 / i+j — 1 i+n—1 i+j i+n—1 

n = ^ (a i+ j + c i+i ) (ofe c fc ) (b k d k ) + (b i+j + d i+j ) JJ(a fe c k ) j [ (6& 4) 



3=0 



k=i 



and 



fc=i+3 - 
i+n—1 



k=i 



k=i+j+l 



i+n—1 



Ki = {di + Cj) I [ (ofc c fc ) + (ft* + dj) ]J[ (6 fc 4) 

fc=j+l fc=i+l 
n—1 / i+j — 1 i+n—1 i+j i+n—1 

+ ^2 [ ( a i+j + c *+i) n (°* n ^ + + n («* ^ n ( & & 4) 

j=l \ fe=i+l k=i+j k=i+l k=i+j+l 

Also define 

Q = - JJ aj Q + JJ 6j dj. 

i i 

Example 4.1. Suppose n = 2. Then 

7i = (a x + c 1 )bidib 2 d 2 + (61 + d 1 )a 1 cib 2 d 2 + (a 2 + c 2 )a x c\b 2 d 2 + (& 2 + d 2 )a x c\a 2 c 2 , 

«i = («i + ci)a 2 c 2 + (&i + d x )b 2 d 2 + (a 2 + c 2 )b 2 d 2 + (6 2 + d 2 )a 2 c 2 . 

Now introduce additional parallel wires with parameters p and — p from 1 to n'. This is 
a special case of the local electrical equivalence for parallel resistors. (Here the parameters 
p and — p should be considered formally, instead of as nonnegative real numbers.) We 
may perform Y — A operations to move the parameter p through the resistor network, as 
shown in Figure [51 










Figure 5. Calculating the electrical R- matrix via certain "virtual" elec- 
trical networks. 



Lemma 4.2. There is a unique non-zero parameter p, which is unchanged after moving 
through one revolution. 

Proof. Let us denote by pi the value of the "extra edge" after 2i star-triangle transfor- 
mations. (The fourth network in Figure [5] shows the location of p\.) We claim that the 
parameter pi is a ratio of two linear functions in the original p = pq. This is easily verified 
by induction using the star-triangle transformation. Thus after 2n star-triangle transfor- 
mations the equation p = p n we obtain is either a linear or a quadratic equation, and zero 
is clearly one of the roots. Therefore there is at most one non-zero solution, and as we 
shall see soon in the proof of Theorem I4.3[ a solution indeed exists. □ 
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Suppose we perform the sequence of Y — A operations of Figure [5j using the parameter 
p of Lemma 14.21 (which we still have to prove exists). Then as illustrated in the final 
diagram of Figure |5j the "extra edges" with parameters p and — p can be removed via 
the local electrical equivalence for parallel resistors. We define the electrical R-matrix to 
be the transformation R(ai,bi,Ci,di) = (a-,^,c^,c^) induced on the edge weights by this 
transformation. 

Theorem 4.3. The electrical R-matrix is given by 

l T i i T i ,/ r i+l ,/ T i+1 

a i = c i = b i = T~ d i = j — • 

Proof. We claim that the parameter p = p Q = Q/kx is the parameter of Lemma 14.21 
Define Pi to be the parameter after i pairs of A — Y and Y — A operations, and qi to be 
the parameter obtained from pi by one A — Y operation. Then 

qi = r 1 - i /Q Pi = Q/ni-i- 
Indeed, calculating by induction 

(ax-i + ci_i)pi + (ax-i cx-i) (ax-i + c X -i)Q + (a X -i cx-i)Kx-i r i-i 



Pi Q Q 

and similarly for pi = bx-idx-i/ \qi-x + + We may then calculate that the 

transformation is given by 

a x _iPi + Cx-iPi + a x _iCi (ax-i + C\-i)Q + a X -iCx-iKx-i Tx-i 



ax-i 0-X-i K l-i CLi—iKx- 



and similarly c\ ,■ = T1 - — . We also have 

,/ qi-xdx-i r 2 -idx-i t 2 - 



qi-x + dx-i + bx-i T 2 -i + (dx-i + h-i)Q bx-iKx-i 
and similarly d'i_i = T2 ~* 1-8 ■ □ 



Proposition 4.4. The electrical R-matrix does not depend on where you attach the extra 
pair of edges, and is an involution. 

Proof If one applies the Y — A transformations to successively push through two edges, 
with weights q and — q negative of each other, then there is no net effect on the weights 
of the other edges involved. Also, the weights on the two edges pushed through remain 
negatives of each other. Thus if one pushes through p to perform the electrical i?-matrix, 
and then a —p, the latter will perform a transformation that will undo the first one. 
Furthermore, a working choice of parameter p at one location this way yields a working 
choice of the parameter at any other location. The uniqueness in Lemma I4T21 implies that 
all resulting electrical i?-matrices are the same. □ 

4.2. Electrical .R-matrix satisfies Yang-Baxter. Let us now consider the network 
iV(3). The procedure of Section H~T1 gives two different electrical i?-matrices acting on 
iV(3): by acting on the part of the network involving vertices M.f\ , or by 

acting on the part of the network involving the vertices , , M- 3 ^ . We denote 
these i?-matrices by R 12 <8> 1 and 1 <g> i? 2 3 respectively. 
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Theorem 4.5. The electrical R-matrix satisfies the Yang-Baxter equation 

(#12 <g> 1) o (1 <g> #23) o (R 12 g> 1) = (1 <g> R 23 ) o (R 12 ® 1) o (1 <g> R 23 ). 

Proof. First we note that to perform the electrical i?-matrix, we can either add extra 
horizontal edges with conductances p and — p between 1 and n', or we could split the 
vertex into three vertices and add vertical edges with conductances p and — p between 
them (see Figure EJ). 



To perform the sequence of .R-matrices (R\2 <8> 1) o (1 (g> #23) ° (-R12 <8> 1), we will add 
horizontal edges for the first and third factor, but add vertical edges for the second 
factor. Let the corresponding weights be p, q, r as shown in Figure [7| According to 



Lemma 14.21 and Theorem 14.31 and Proposition I4.4[ these weights exist and are unique. 
Now, apply the A — Y transformation as shown in Figure |7] to obtain weights p f , q' 
and r', and their negatives on the other side. We claim that if these new weights are 
pushed around the cylinder, they come out the same at the other end. This follows from 
Proposition 12.31 Now if we push the weights p, q and r through and apply the A — Y 
transformations, we obviously get again the weights p', q' and r'. Thus by Lemma 14.21 
and Theorem 14.31 and Proposition I4.4[ while pushing p', q' and r' through we are applying 
(1 <g> #23) (-R12 <E> 1) ° (1 <g> #23)- Since, by the electrical tetrahedron relation (Proposition 
12.31 ) the two results are the same, the claim of the theorem follows. □ 




Figure 6. 




Figure 7. 
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Corollary 4.6. The electrical R-matrix gives an action of the symmetric group S m on 
N(m). 

Proof. Follows from Theorem 14.51 and Proposition 14.41 □ 

5. Electrical ASW factorization 

5.1. Universal response matrix. Let G be a cylindrical electrical network. Thus G 
is an electrical network embedded into a cylinder so that the intersection of G with the 
boundary of the cylinder is exactly the boundary vertices of G. Let C7(oo) denote the 
universal cover of G. It is an infinite periodic network embedded into an infinite strip. 
Given finite sets V\ C Vi C . . . of vertices of G(oo) which eventually cover all vertices 
of G(oc), we obtain a sequence of truncations G(l), (7(2),... of G(oo) as follows. We 
let G(N) be the subgraph of G(oo) consisting of all edges incident to a vertex in Vjv- 
Furthermore, we declare a vertex of G(N) internal if it lies in Vjy and is internal in G(oo). 
Each G(N) is a finite planar electrical network. 

We suppose that the boundary vertices of G(oo) are numbered Z = {. . . , —2, —1, 0, 1, 2, . . .} 
on one side of the boundary, and by U — {. . . , —2', —1', 0', 1', 2', . . .} on the other side 
of the boundary. (We will always picture the infinite strip as vertical, with vertex labels 
increasing as we go downwards.) We define the universal response matrix of G to be given 
by Cij where 

C^j = lim LijiGiN)) 

where i,j denote vertices of G(oo). For sufficiently large N, any two fixed vertices of 
C7(oo) will lie in G(N). These limits exist and are finite, due to following lemma. 

If H is an electrical network, and V a subset of its vertices, the response matrix of H/V 
is the response matrix obtained by declaring all the vertices in V to be interior. 

Lemma 5.1. Assume V C V are two subsets of vertices of an electrical network H, and 
assume i and j are two vertices not contained in V . Then L'^ in the response matrix of 
H/V is at least as large as Lij in the response matrix of H/V. 

Proof. The Schur complement with respect to some set V can be taken as a sequence of 
Schur complements with respect to single vertices in V in some order (see for example 
\CM\ (3.7)]). Thus, it would suffice to prove the statement for V — V consisting of a 
single vertex. In this case the claim is obvious however, since off-diagonal entries of a 
response matrix are nonnegative, and diagonal entries are non-positive. □ 

Theorem 5.2. The matrix Cij is a well-defined infinite periodic matrix, which does not 
depend on which truncations G(N) are taken. 

Proof. There are several parts to this statement. 

The limits Cy exist. Note that for sufficiently large N, Lij(G(N)) can be calculated by 
taking of the network G(oo)jv, which is obtained from G(oo) by declaring that only 
the internal vertices of G(N) are internal in G(oo)n- The network G(oo)n is obtained 
from G(N) by adding extra boundary vertices attached only to boundary vertices of G(N) 
(and by assuming N is large enough, we may assume that these extra vertices are not 
incident to i or j). But Lij(G(N)) is by definition calculated by measuring the current 
flowing through j when vertex i is set to one volt and all other boundary vertices are 
set to zero volts. Since current does not flow between zero volt vertices it follows that 
Lij(G(N)) = Lij(G(oo))]y. By Lemma I5TT1 the sequence Lij(G(N)) N is non- decreasing as 
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N — > oo, since each network is obtained from the previous one by declaring some extra 
vertices internal and taking the corresponding Schur complement. 

The limits dj do not depend on the sequence of truncations. Assume we have two 

different sequences V\ C V 2 C . . . and V{ C V 2 ' C Since we know that each eventually 

covers all vertices, we know that for each i there is a j such that V{ C V- and V- C Vj. 
Then applying Lemma I5TT1 we conclude that the two limits bound each other from above, 
and thus are equal. 

The limits dj are periodic: dj = C(i+ n )(j+ n ). Indeed, take two sequences V\ C V 2 C . . . 
and V[ C V 2 C . . ., one obtained from the other by a shift on the universal cover by the 
period n. We know that they give the same value of dj by the previous part. On the 
other hand, it is clear that the value one of them gives for dj is the value the other gives 
for C(i +n )(j +n ). 

The limits dj are finite. Any truncation of G(oo) can also be viewed as a truncation 
of a finite cover G[m] (obtained by lifting G to a m-fold cover of the cylinder) for a 
large enough m. By Lemma 15.11 the conductance Lij(G(N)) is bounded from above 
by the same conductance in G[m], which in turn is bounded from above by the same 
conductance in the original network G. Indeed, if ji, ■ ■ ■ ,j m are vertices in G[m] that 
cover j, then Lij(G) = J2T=i Lij k (G[m]). This can be seen as follows. Using the linearity 
of the response and periodicity, the sum Yuk=i Lij k (G[m]) measures the current through 
je (for any £) when all vertices if. that cover % have potential 1, and all other vertices have 
potential 0. But projecting onto G by identifying all covers of a vertex we see that this 
current flow is exactly L{j(G). In particular, Ly(G(jV)) < Ly fc (G[m]) < Lij(G) for any 
k. This shows that dj is bounded from above by Lij(G), and thus if the latter is finite, 
so is the former. □ 

Thus for fixed we can approximate dj arbitrarily well by calculating Lij(G(N)) 
for some large N. 

Proposition 5.3. The universal response matrix dj i g invariant under the local electrical 
equivalences of N. 

Proof. For any local electrical transformation in G one can choose a sequence of trunca- 
tions of G(oo) containing completely several occurrences of this transformation. Since the 
conductances in these truncations do not change, the limit is also invariant. □ 

Corollary 5.4. The electrical R-matrix preserves the universal response matrix. 

Proof. The only comment one needs to make is that the entries of the universal response 
matrix are limits of rational functions in the edge weights, and that these rational functions 
make sense formally even when negative conductances are used (as in the derivation of 
the electrical R- matrix). □ 

5.2. Circular and cylindrical total nonnegativity. Let I — {ii < i 2 < ... < 4}, J = 

{ji < 32 < ■ ■ ■ < jk} C [n] be subsets of the same cardinality. Then (/, J) is a circular 
pair if a cyclic permutation of ii,i 2 , ■ ■ ■ ,ik,jk, ■ ■ ■ ,ji is in order. A n x n matrix M is 
circular totally-nonnegative if the minor det(M/ ) j) is nonnegative for every circular pair 
(J, J). Curtis, Ingerman, and Morrow [CIMJ show that the response matrices of circular 
planar electrical networks are exactly the set of circular totally-nonnegative symmetric 
matrices for which every row sums to 0. (Note that the response matrices in [CIMj are 
the negative of ours.) 
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Let us extend this to cylindrical electrical networks. Put the total order • • • — 1 < 
< 1 < 2 < • • • < 2' < V < 0' < -1' < • • • on Z U II. Let J, J C Z U II be 
two ordered subsets of the same finite cardinality. Then (J, J) is a cylindrical 'pair if 
a cyclic permutation of i 1; i 2 ? • • • ^kyjk, ■ ■ ■ ,ji is in order. A matrix M with rows and 
columns labeled (and ordered) with ZUZ' is cylindrically totally-nonnegativeii the minor 
det(M/ j) is nonnegative for every cylindrical pair (/, J). 

Proposition 5.5. Suppose C = (Cij) is the universal response matrix of a finite cylin- 
drical electrical network. Then C is cylindrically totally-nonnegative. 

Proof. For each fixed cylindrical pair (/, J), and sufficiently large N, the truncation G(N) 
is a finite circular electrical network including all the boundary vertices in I and J. 
But then det(£/ i j) = lim N ^. QO det(L I; j(G(N))) > 0, using Curtis-Ingerman-Morrow's 
result. □ 

Remark 5.1. The converse to Proposition I5.5[ namely, which cylindrically totally non- 
negative matrices are realizable as universal response matrices, is more subtle. Let 
. . . , j-i,jo, ji, j2, ••• be the lifts to the universal cover of a particular vertex in a finite cylin- 
drical electrical network G. Then for a fixed i, the (doubly-infinite) sequence = Ci t j k 
must satisfy certain recursions or convergence properties. In the different but closely re- 
lated setting of total nonnegative points of loop groups, the correct property is to ask for 
the generating function of to be a rational function (see | LP10+j Theorem 8.10]). 

5.3. ASW factorization for networks N(m). We now assume we are given a network 
G = N(m). The vertices have been labeled so that if we take the "low" edge (from M( 
to Mf +1) ) at every step starting from i we will end up at i' . Note that a shortest path 
from one side of the cylinder to the other consists of exactly m edges. 

From now on we consider groves in the universal cover G(oo) of G (or in the truncations 
G(N)). The boundary partition of such a grove would be a planar partition of Z U Z' 
arranged on the two edges of an infinite strip (or in the truncations of this) . 

Lemma 5.6. Suppose i and j' are can be connected by a path with m edges. There exists 
an integer M such that there are no groves V with the properties 

(1) there is a grove component ^{ij>} with boundary vertices {i, j'} and which uses an 
edge below (resp. above) any of the shortest paths from i to j' . 

(2) there are grove components with boundary vertices {i + 1, (j + 1)'}, . . . , {i + M, (j + 
M)'} (resp. {i - 1, (j - 1)'}, . . . ,{i — M,(j — M)'}). 

Proof. For a grove component ^{i t j>} let us call bad the edges it contains that are below 
the lowest path from i to j'. Assume the grove component ^{ij'} has bad edges, and 
furthermore without loss of generality assume that it has a bad high edge. The case of a 
bad low edge is similar with the left and right sides of the network swapped. Assume k 

(k) 

is the first index such that bad high edge has one of the M- as its right endpoint. We 
claim that the F^+i^+iyy component has a bad high edge with right end having index 
k — 2 or smaller. Indeed, consider the unique path from i + 1 to (j + 1)' inside T^i^iyy. 
In order to avoid touching the bad edge of ^{i,f} it has to turn, diverting from the lowest 
shortest path from i + 1 to (j + 1)'. The first time it thus diverts gives a desired high 
edge. 
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Figure 8. A bad high edge in the blue grove component forces an earlier 
bad high edge in the red grove component. 

Now, the index of the first bad high edge cannot decrease indefinitely, and in fact one 
sees that the statement of the lemma holds for M > m/2. The case of edges above the 
highest shortest path is similar. □ 

Define the radii Rk for k — 1,2, ... ,m by 

rik — T^f 

where we have denoted the conductance of the edge joining two vertices v , v' by w vv i. 

Lemma 5.7. Let (Pj, P 2 , . . . , P n ) be an n-tuple of consecutive shortest paths, where Pi 
connects i + k to (i + k — 1)' for each i, and some fixed k. Suppose that 

(1) all the Pi have the same shape; thus they high or low edges respectively at the same 
points along the path. 

(2) R 1 > R 2 > ••• > R m . 

Then the total weight wt(Pi U P 2 U • ■ ■ U P n ) is maximized exactly when the highest path 
is taken. 

The weight of a subgraph is simply the product of its edge weights. 

Proof. All the shortest paths are connected by switching from one side of a rhombus to 
the other, see Figure [91 If this parallelogram involves vertices with upper index k — 1, k 
and k + 1, then the higher path has greater weight exactly when R^ > Rk+i- D 

Let Tk be the partition of Z U U with parts of size two {k, (k — 1)'} for all k 6 
{2, 3, ... , K}, parts of size two {k, k'} for k £ {— K, —K + 1, . . . , 0}, and all other parts 
are singletons. Thus in particular, 1 is a singleton. We shall denote by o~k the partition 
of Z U Z' obtained from by placing 1 in the same part as {0, 0'} to get a single part 
{1,0, 0'} of size three. 

We shall suppose that the truncation G(N) includes the boundary vertices — N, —N + 
1, . . . , 1, 2, . . . , iV and also the boundary vertices —N', (-N + 1)', . . . , 1, 2, . . . , N'. For 
N 3> K the partitions tk and o~k naturally gives rise to partitions of the boundary 
vertices of G(N), where boundary vertices of G(N) which are not boundary vertices of 
G(oc) are all considered singletons. We will still denote these partitions of boundary 
vertices of G(N) by tk and o K . 
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Figure 9. Here n = 2; which of the two groves has larger weight depends 
on which of the two radii ac/bd and fh/eg is bigger. 

Theorem 5.8. Suppose that Ri > R2 > • • • > R m - Then 



where a is the weight of the high edge connected to the vertex 1 of G. 

Proof. Let v denote the vertex connected to both and 1. Let e be the edge joining 1 to 
v, so that e has weight a. Let e' denote the other edge incident to 1. 

To approximate the LHS, we shall assume we have chosen N ^> K ^> M. Let T be 
a grove with boundary partition either o~k or tk- By Lemma 15.61 the grove component 
Fk-m,(k-m-i)i cannot extend either above or below the set of edges contained in shortest 
paths from K — M to (K — M — 1)'. In particular, there is a bound on the number of 
choices of Tk-m,{k-m-i)'i n °t depending on K or N. In the following, we shall assume 
that we have fixed such a choice for T k-m,{k-m-\)' ■ 

Let T be a grove with boundary partition tk- By Lemma [5761 the grove component r 0i o' 
cannot extend above the (unique) shortest path from to 0'. In particular, the vertex v 
must lie in r 0) o'- Thus the edge e is never used in T, and T U {e} is a grove with boundary 
partition ax- Also by Lemma I5"75| the grove component T-m,(-m)> cannot extend either 
above or below the (unique) shortest path from — M to (— M)' , and therefore must be 
exactly this shortest path. Similar observations hold for a grove with boundary partition 
ok- In particular, the part of the grove above F_ M / M y and the part below are essentially 
independent. 

Let us denote by A T (resp. A a ) the set of groves with boundary partition tk (resp. ctk) 
and by wt(A T ) (resp. wt(y4 CT )) the total weight of that set of groves. 

Let T be a grove with boundary partition ok- First we observe that all but a constant 
number of grove components I\ ($_]V are shortest paths. Furthermore, as % goes from 2 to 
K, the shapes of the shortest paths are locally constant, and can only change when we 
encounter a grove component which is not a shortest path. Furthermore, as we go down, 
the shape of the shortest path can only become lower. We shall call this part of T the 
lower half of the grove. Note that there are m different shapes Si, S2, ■ ■ ■ , S m of shortest 
paths, listed from highest to lowest, see Figure [101 Let A' a C A a denote the subset of 
groves where r 10 ,o' uses the edge e, and let A" denote the subset of groves where r 10 ,o' 
does not use the edge e'. Given e > 0, we shall now show that for sufficiently large K, 
one has wt(A")/wt(^) < e. 




= a 



16 



THOMAS LAM AND PAVLO PYLYAVSKYY 




Figure 10. The shapes of shortest paths connecting i + 1 to i' . 



Let r G A". Then the shortest paths which occur in the lower half of T can only use the 
shapes S2, S3, ... , S m . Pick K sufficiently large that we are guaranteed to have Cnm 2 /e 
grove components in the lower half which are shortest paths, where C is a constant (not 
depending on K or N) we shall describe below. Then there is some Si (say pick the 
least such i) which occurs at least Cnm/e times. Define a set of new groves {T'j \ j = 
1,2,..., m/e} by: 

(1) removing jn of the Si shaped paths, and shifting the lower half of T downwards 
to fill in the removed area; 

(2) replacing the grove components Y_ M ^_ M y, r_ A f+i,(-A/+i)'? • • • ? rL^-iy with short- 
est paths; 

(3) replacing r ,o',i with the shortest path from to 0' union the edge e; 

(4) adding jn — 1 new shortest paths Y' 2V , . . . , Y'- n u n ^iy which are Si-shaped; 

(5) adding one extra "transition" grove component Y'- n+1 ^ n y which is ^-shaped 
above, but which correctly fits with r 2) i/ below. 

Figure [Til illustrates the procedure. In this case m — 9, n — 2, j — 1, i — 7. The Ti^o' 
component of the groves is shown in red, the part shifted down is shown in green. The 
two removed groves of shape Si and the two new groves of shape Si that replaced them 
are shown in brown. 

Note that T'j belongs to A' a . Since apart from the shortest paths the rest of the grove 
has only been modified at a bounded number of edges, by Lemma 15.71 there is a constant 
C\ not depending on K,N, or e such that wt(r^) > Ciwt(r). Furthermore, each grove 
in A' a can occur in this way in at most C^vn different ways. The C2 counts the possible 
grove components Y_ M r M y, r^+i^-iw+i)') • • • > T-i^-iy and r 0j0 ',i' of which there is a 
universal bound for. The m ways count the possible i such that the shortest paths Si are 
being replaced by S±. Setting C = G\G%, we conclude that wt(A^) < ewt(A' a .). For every 
G A T , the grove U {e} lies in A a , and this map is an injection from A T to A a which 
changes the weight of each grove by exactly a. On the other hand every r G A' a is in the 
image of this map. It follows that a < wt(A a )/wt(A T ) < a(l + e). Letting K,N — >■ 00, 
we obtain the statement of the theorem. □ 

The partition tk is a planar partition of ZUZ' of the type described in Theorem 13. 1( 2). 
where all but finitely many vertices are isolated. Define Vr(r) by taking the polynomial 
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Figure 11. Creation of a new grove 1^ from T. 

of Theorem 13.1( 2) and replacing Lij(G) by Cij. Then we have 

Vr(r) = lim Pt (t)g(n) 

and in particular, Vr[f) can be approximated arbitrarily well on some G(N) for very 
large N. 

Unfortunately, this is not the case for the partition cr^, which contains a part of size 
three, for which Theorem 13.1( 2) cannot be applied. Instead one has 

Pr(a K ) G(N) =p N (L tJ (G(N))) 

for the sequence of polynomials pn of Theorem 13 .1( 1). These polynomials depend on <jk, 
which is suppressed from the notation. 

Proposition 5.9. Suppose N > N' . Then p^ ~Pn' is a polynomial in the ; s such that 
every monomial involves some Lij where N > \i\ > N' and j G [—K,K]. Furthermore, 
there is some constant ck such that the coefficient of each monomial in pn is less than 
c K . 

Proof. The first statement follows from Lemma 13.31 and our choice of The second 
statement is Lemma [3.41 □ 

The assymmetry of the roles of % and j in Proposition 15.91 is accounted for by the fact 
that is symmetric. That is, we treat Ly = Lji as the same variable. Note that we 
already know liniTv^oo P t ( t k)g{n) approaches a limit, and it follows from Theorem 15.81 
that liniAr_ s . 0O Py(o~k)g(n) approaches a limit as well. 

Lemma 5.10. For each K and e > ; one can find some A such that 

| lim Yi(a K ) G{N) - p A (C t j)\ < e. 
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Proof. Fix K. We first show that there is A such that for all N > A we have 

\ PN {L l0 {G{N))) - p A {L %3 {G{N)))\ < e/3. 

It is known that the polynomials p^ have degree IK + 1, see Section [3) By Proposition 
15.9} every monomial in p^ — p A has a factor where |i| > A and \j\ < K, and has 
coefficient < ck- Thus 

/ 

| PJ v(^(G(iV)))-p A (^(G(iV)))|< 



Eiij (max(L a6 )) 2X c#. 
I nun 

\M>A li|<x 



By the proof of Theorem 15.21 we know that for each i we have . A? < 00 • Thus it is 
possible to find A large enough that 




< 



3(max a , b (L ab )) 2K c K 



giving us 

(2) \p N (L t3 (G(N))) - p A (L tJ (G(N)))\ < e/3. 

Since only finitely many Ly's appear in p^, for sufficiently large iV we have 

(3) VAiC^-pAiL^GiN))) <e/3. 
Furthermore, for sufficiently large iV we have 

(4) | Pr (cr^) G(Ar) - lim Yr{a K ) G ( N n\ < e/3. 

Finally, we combine the three estimates (E]),(j3]), @ and use Pt(o-k)g(n) = PN{Lij(G(N))). 

□ 

Theorem 5.11. FixG, satisfying Ri > i? 2 > • • • > -R m and fix one of the edges connected 
to one of the vertices . . . , —2, — 1, 0, 1, 2, . . .. There is a sequence of polynomials (depending 
on the universal response matrix of G), p Al ,p A2 , ■ ■ ., and qi, q 2 , ■ ■ ■ such that 



lim 



K->oo q K {C>ij) 
where a is the weight of the chosen edge. 

Proof. By symmetry, it is enough to establish the formula for the high edge connected to 
1. The polynomial qx is the one associated to from Theorem 13.1( 2). The polynomial 
Pa k is chosen via Lemma [5. 101 so that \p Al< (Cij) — Ymi n ^ 00 ~Pi:(o-k)g{n)\ < £k, where ex 
is chosen so that q ^..^ ~ ► as K — )■ oo. Finally, we apply Theorem 15.81 □ 

6. The Inverse Dirichlet-to-Neumann problem on a cylinder 
6.1. Solution to inverse problem for the networks N(m). 

Lemma 6.1. Assume a cylindrical network Y = N(m) is obtained by concatenating two 
networks X = N(l) and X' = N(m — 1). Then knowing X and the universal response 
matrix ofY, one can recover the universal response matrix of X' . 



INVERSE PROBLEM IN CYLINDRICAL ELECTRICAL NETWORKS 



19 



Proof. Assume the conductances in X are for the high edges and 6j for the low edges. 
Concatenate Y with a network Z = N(l) with (virtual) conductances — Oj for low edges 
and —hi for high edges. We claim that the response of the resulting network is equal to 
that of X'. Indeed, connect the opposite vertices of Z and X by edges with conductances 
L and — L, without changing the response. Changing resulting triangles into stars using 
the Y — A transformation, and letting L — > oo, we see that there is an infinite conductance 
between opposite vertices and zero conductance between other pairs. Thus the two iV(l) 
networks effectively cancel each other out, and we are left with a network with the same 
response as X' . See Figure [12j □ 




Figure 12. Attaching Z to Y effectively removes the leftmost layer X. 

Theorem 6.2. There are generically ml sets of edge conductances which produce a given 
universal response matrix for the network N(m). All the solutions are connected by the 
S m -action via the electrical R-matrix. 

Proof. It is easy to see from Theorem 14.31 that the electrical .R-matrix swaps the radii 
Rk- Assume we have a solution for conductances in N(m) with given universal response 
matrix. Apply the electrical .R-matrix to reorder the radii Rk in non-increasing order. 
Then Theorem 15.111 allows us to recover the conductances in the leftmost N(l) part of 
the network. In particular, these conductances are the same for any solution. Once we 
know that, we can use Lemma lBTTl to recover the universal response matrix of the remaining 
N(m — 1) part of the network. Then we repeat the procedure. We see that once we require 
the radii to form a non-increasing sequence, the conductances are recovered uniquely by 
Theorem 15. Ill Therefore all other solutions can be obtained from that one by action of 
S m , which is what we want. In the generic case when all radii are distinct, the orbit has 
size ml. □ 

Remark 6.1. In the language of | LP10+ |, Theorem 16.21 says that the group S m generated 
by electrical .R-matrices is exactly the monodromy group of the network N(m). It is the 
group acting on the edge weights of N(m) obtained by transforming the network via local 
electrical equivalences back to itself. 

6.2. Conjectural solution to general case. It is convenient to describe the general 
answer we expect using the language of medial graphs, see for example [CIMl IdVG Vj . 
Draw wires through the electrical network so that they pass through each edge and connect 
inside each face as shown in the first two pictures in Figure [121 The third picture shows an 
example of an electrical network and its medial graph. Note that the medial graph always 
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has four-valent vertices, and that the wires of the medial graph "go straight through" 
each vertex. 



Figure 13. Constructing the medial graph: the black edges belong to the 
electrical network, and the blue lines constitute the medial graph. 

The medial graph of the networks N(m) looks like 2n horizontal wires crossed by m 
cycles, as shown in Figure HU 



In the terminology of |CIM] a circular planar electrical network is critical if its medial 
graph avoids lenses, which is equivalent to saying that every pair of wires crosses as few 
times as possible, given their respective homotopy types. If G is a cylindrical electrical 
network, we say that G is critical if the universal cover of G satisfies this condition; 
namely, the medial graph has wires which cross as few times as possible. 

Let us call a cylindrical network canonical if the medial graph of its universal cover has 
the following form. First, there are three kinds of wires: (I) some wires connect points 
on opposite boundaries, (II) some wires connect points on the same boundary, and (III) 
some wires do not intersect the boundary at all and correspond to (simple) cycles around 
the cylinder. Secondly, we require that the third kind of wires do not intersect the second 
kind, and furthermore, all points of intersection of wires of the first kind with themselves 
happen strictly before they intersect wires of the third kind. An illustration is given in 
Figure [T5j 

Conjecture 6.3. 

(1') Any cylindrical electrical network can be transformed using local electrical transfor- 
mations (those in Section \2J% and the electrical R-matrix) into a critical cylindrical 
electrical network. 

(2') Any two cylindrical electrical networks G and G' with the same universal response 
matrices are connected by local electrical equivalences. Furthermore, if G and G' 
are both critical, then only star-triangle transformations, and electrical R-matrix 
transformations are needed. 





Figure 14. (N{m) and its medial graph. 
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Figure 15. The (universal cover of the) medial graph of a canonical network. 

(3') If a cylindrical electrical network is critical canonical, then the conductances cor- 
responding to all crossings involving wires of types (I) and (II) can be recovered 
uniquely. The conductances corresponding to crossings of wires of type (I) and 
type (III) can be recovered up to the electrical R-matrix action. 

(5') The space X of universal response matrices of cylindrical electrical networks has 
an infinite stratification by X = UCi where each Ci ~ M> x K> is a semi-closed 
cell that can be obtained as the set of universal response matrices for a fixed critical 
network with varying edge weights. 

For the (missing) cylindrical analogue of (4) of the Theorem in the introduction see 
Section 15.21 

Another way to phrase Conjecture 16.3( 3') is that the monodromy group of a critical 
canonical cylindrical network is a symmetric group, generated by electrical i?-matrices. 
See Remark [6.11 

Let us explain the semi-closed cells in Conjecture 16.3( 5') . Let G be a critical canonical 
cylindrical electrical network. Some edge weights can be recovered uniquely and these each 
give a M>o in the parametrization. The remaining part of the network is essentially one of 
the networks N(m), whose edge weights can be recovered uniquely up to the electrical R- 
matrix action (Theorem l6.2p . So the response matrices would be parametrized by the orbit 
space (IR'^q )/S m . However, we can pick a distinguished element in each orbit: namely the 
one where the radii Rk are non-increasing. The corresponding response matrices would 
then be parametrized by R m G M >0 , R1 — R2, R2 — R3, ■ ■ ■ , R m -\—R m G R>o together with 
some collection of edge weights which can be freely chosen in M>o- Thus the universal 
response matrices of critical canonical cylindrical electrical network is parametrized by a 
semi-closed cell C ~ M> x M> . 
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